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Global fixed point proof of time-dependent density-functional the- 
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Abstract. - We reformulate and generalize the uniqueness and existence proofs of time-dependent 
density-functional theory. The central idea is to restate the fundamental one-to-one correspon- 
dence between densities and potentials as a global fixed point question for potentials on a given 
time-interval. We show that the unique fixed point, i.e. the unique potential generating a given 
density, is reached as the limiting point of an iterative procedure. The one-to-one correspondence 
between densities and potentials is a straightforward result provided that the response function 
of the divergence of the internal forces is bounded. The existence, i.e. the v-representability of 
a density, can be proven as well provided that the operator norms of the response functions of 
the members of the iterative sequence of potentials have an upper bound. The densities under 
consideration have second time-derivatives that are required to satisfy a condition slightly weaker 
than being square-integrable. This approach avoids the usual restrictions of Taylor-expandability 
in time of the uniqueness theorem by Runge and Gross [Phys. Rev. Lett. 52, 997 (1984)] and of 
the existence theorem by van Leeuwen [Phys. Rev. Lett. 82, 3863 (1999)]. Owing to its general- 
ity, the proof not only answers basic questions in density-functional theory but also has potential 
implications in other fields of physics. 
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Time-dependent density-functional theory (TDDFT) 
provides an exact reformulation of the time-dependent 
quantum many-body problem in terms of an effective non- 
interacting problem that is functionally dependent on the 
one-particle density. Since the effective problem is nu- 
merically much more tractable than the original one the 
method has found widespread use in different fields of 
physics and chemistry [IH3] • The cornerstone of TDDFT 
' is the existence of a one-to-one mapping between densi- 
ties and external potentials for a given initial state. The 
standard proofs that establish this mapping by Runge and 
Gross [4] and van Leeuwen [5] restrict the set of allowed po- 
tentials to those that have a Taylor-expansion around the 
initial time. In this work we extend these proofs to allow 
for general finite potentials without restrictions on their 
time-dependence. Hence we can, for example, include po- 
tentials of the temporal form \ft or exp (— 1/t) around 
t = 0. This was until now only proven to be allowed for 
the linear response case [2], dipole fields [7] or lattice sys- 
tems [§]. The proof of existence [5] also requires Taylor- 



expandability of the density. This issue is discussed for the 
one-dimensional case in [B] . In the present formulation we 
do not need this property. The second time-derivative of 
the density is only required to satisfy a condition which is 
slightly weaker than being square-integrable. The central 
idea to go beyond the Taylor-expandability restrictions 
is to reformulate the fundamental one-to-one correspon- 
dence between densities and potentials as a global fixed 
point question for potentials on a given time-interval. We 
prove that the fixed point, i.e. the potential that gener- 
ates a given density, is unique and exists under certain 
assumptions. 

Let us outline the general problem. We consider a many- 
body Hamiltonian in atomic units of the form 



H([v],t)=T + V([v],t) + W 



(1) 



where the operator T = — $^ i=1 is the kinetic term, 
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W = | J2ijtj w(\ri — r_y|) is the two-particle interaction, 

N 

V([v],t) = ^2v(rit) = / drn(r)v(rt) (2) 

i=i 

is the external potential and 

JV 

n{v) = Y J 5{v-v l ) (3) 

i=l 

is the density operator. The interaction potential w(\r — 
r'|) can be any locally square- integrable potential with 
a Coulomb singularity or less singular, though it is usu- 
ally assumed to be purely Coulombic. For a given initial 
state l^o) at time to the solution of the time-dependent 
Schrodinger equation (TDSE) 

id t M[v],t))=H([v],t)\*([v},t)) (4) 

for a potential v (including gauge fixing) provides a map- 
ping from potentials to wave- functions v t— > \^([v],t)). 
The wave-functions are as usually assumed to be square 
integrable and in the self-adjoint domain of the Laplacian, 
i.e. V 2 V> is also square integrable. Due to the mapping 
of potentials onto wave functions the expectation value 
of any operator O becomes a functional of the external 
potential, i.e. 

0([v],t) = (*([v],t)\6m[v],t)). (5) 

In particular the density n([v],rt) is a functional of the 
potential and the initial state. The central statement on 
which TDDFT is based is that this mapping v H> n is one- 
to-one and thus invertible, i.e. for a given density there is 
only one external potential generating this density by solu- 
tion of the TDSE for a given initial state l^o) • This implies 
that the wave-function and all observables can be consid- 
ered as functionals of the density instead of the poten- 
tial. The statement is independent of the initial state and 
the two-particle interaction. This opens up the possibility 
that the density of an interacting system can be repro- 
duced by a unique effective potential in a non-interacting 
system. This forms the basis of the Kohn-Sham scheme 
that makes TDDFT computationally practical. Before we 
outline the general proof we reinvestigate the following 
equation, which is of fundamental importance to TDDFT 

m 

- V- [n([v],rt)Vv(rt)] = q([v], rt) - <9 t 2 n(M, rt). (6) 

It can be obtained by applying the Heisenberg equation of 
motion for the density operator twice. In this equation 

q([v],rt) = (*(H,i)|ft(^f w (r) + Wi(r))|*([«],t)> (7) 

where summing over multiple indices is implied and 

, N 



- i 9 fe a ; n(r) (8) 

N N 

m(r) = ^(r-r^XjMr-r,-!)) (9) 

»=1 j¥=i 

represent the momentum-stress tensor and internal force 
density of the system [9]. Here d k,i and a Yi are par- 
tial derivatives of particle i that act to the left and right 
respectively. If we replace n([v],rt) in eq.© by a given 
density n subject to the following conditions on the initial 
density and its time-derivative at time to 

n(rt ) = (# |ft(r)|*o) (10) 
d t n(rt ) = -<*o|V-j(r)|* > (11) 

where j(r) is the usual current-density operator with 

Jk(r) = ^ fe S ( r ~ ~ ^ k ^ r ~ r *)) ' ( 12 ) 

then eq.© becomes 

- V • [n(rt)Vu(rt)] = q([v],rt) - d?n(rt). (13) 

This is a nonlinear equation for v which needs to be solved 
with specified boundary conditions (this amounts to fixing 
a gauge for v). If we propagate the TDSE with initial state 
l^o) and with a potential v that is a solution to ea. dTHl) 
then for this potential clearly also the local force equation 
([6]) will be satisfied with the same initial conditions (fl0|) 
and (fTT|) . Subtracting eq.(fl~3]l from ((6]) then yields the 
equation 

dtP(rt) - V • [p{rt)Vv(rt)} = (14) 

for the density difference p(rt) = n([v],rt) — n(rt) with ini- 
tial conditions p(rto) = dtp(rto) = 0. Note, that the solu- 
tion has further to fulfill J drp(rt) = 0. The unique solu- 
tion of eq. (fT4|) with these initial conditions is p(rt) = and 
hence n(rt) = n([v], rt), i.e. the density in ea. p^)) is iden- 
tical to the one that is obtained from time-propagation of 
the TDSE with the solution of eq. (fl~3"|) . By making differ- 
ent choices for the density n(rt) in eq. (fT3)) we can deduce 
some important consequences of this result. If we choose 
n(rt) — n([u],rt) to be the density obtained by time- 
propagation of the TDSE with potential u and the same 
initial state j^o) then we must have n([u],rt) = n([v],rt) 
where v is a solution of eq. (|13|) . The uniqueness of a poten- 
tial v for a given density n (the Runge-Gross theorem) is 
thus equivalent to the uniqueness of the solution of eq. (fH{|) , 
i.e. u = v. If we choose n(rt) to be the density obtained 
by solving a TDSE for a system with different two-particle 
interactions W' and with different initial state |$o) then 
the existence of a solution to eq. (fT3|) implies that the same 
density can be reproduced by a potential v in our system 
with interaction W and initial state \^o), i-e. it is v- 
representable in our system. For the special case W = 
this amounts to reproducing the density of an interacting 
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system within a non-interacting system, which is known 
as the Kohn-Sham construction and forms the basis of 
virtually all applications of TDDFT. The key question, 
which is crucial for the whole foundation of TDDFT, is 
thus whether a solution to eq. fTB")) is unique and exists. 
Existence and uniqueness have indeed been established 
4,5! under the restriction that the potential v is Taylor- 
expandable around the initial time. The main goal of this 
paper is to present a proof that lifts this restriction. 

Before going into details we first give the main idea. 
For a potential vo(rt) and initial state |^o) at time to we 
can propagate the TDSE in a given time interval [to,T] 
and construct the function q([vo],rt) from eq.©. This 
provides a mapping 

V : v k> q[v ] (15) 

(see figH]). Let us fix a density n(vt) subject to the condi- 
tions (fTOf and (fTTj) in terms of the initial state |\&o)- Then 
with the inhomogeneity q[vo] — d 2 n we can solve eq. (|13[) 
for a new potential V\, i.e. 

- V • [n(rt)V«i(rt)] = q([v ],rt) - d 2 t n{vt) (16) 

with given boundary conditions (this amounts to fixing 
a gauge for V\). We point out that this is now a lin- 
ear Sturm-Liouville problem parametrically dependent on 
time and as such not an evolution equation. This provides 
us with a second map 

V : q[v ] ^ Vl (17) 
(see figH}. The combined map 

F[v ] = (yoV)[v ]=v 1 (18) 

maps our original potential vo to a new one v\. If for 
some potential v we have v — T[v\, i.e. v is a fixed point 
of the mapping J 7 , then we satisfy eq. (|13[) . Consequently, 
the question whether a solution to eq.(Q~3|) exists and is 
unique is equivalent to the question whether the mapping 
T has a unique fixed point. This is exactly what we will 
show in this work. Our proof is based on the following 
inequality 

- J"[«o]||a < - «o||a (19) 

with a < 1 and where || • || Q is a norm dependent on a pos- 
itive parameter a on the space of potentials. Comparable 
norms are commonly used in the solution of initial- value 
problems [10J . This inequality can be derived in two steps. 
In the first step we prove that if two potentials are close 
in norm then also the observables 0([u], t) calculated from 
them are close. In particular for O = q(r) we prove 

C 

\\q[vx] - q[v ]\\ a < - v \\ a (20) 

Va 

with C a positive constant. In the second step we prove 
that 

II^M - J>o]|U < D\\q[ Vl ] - q[v ]\\ a (21) 




Fig. 1: The potential-potential mapping T of eq. (|18|) as com- 
position of the mappings V and V. 



for a positive constant D. The combination of these two 
statements then yields eq. (fT9|l where a — CD/ y/a and we 
can choose y/a > CD. It remains to prove Eqs. ([20|) and 

(ZD- 

Before we do so, we point out conditions on the set of 
densities and potentials and some details of the mathemat- 
ics involved. All derivatives and solutions of the differen- 
tial equations are defined in the weak sense, i.e. distribu- 
tionally |10j . In order for the Sturm-Liouville boundary 
value problem of mapping V to be mathematically well 
defined we have to restrict our considerations on a finite 
region f2 C R 3 with boundary dfl. This is not a serious 
restriction since this region can be arbitrarily large, e.g. as 
big as the laboratory or the solar system, and we can ap- 
proach the continuum limit to any accuracy that we desire. 
Note that this restriction is also needed in the original ex- 
istence proof [5] as was discussed in [TT]. The Runge- Gross 
theorem [J] also uses the Sturm-Liouville equation. As was 
already pointed out in [12] and discussed to some extent 
in |13j the proof needs further asymptotic conditions on 
the densities and potentials. These asymptotic require- 
ments are most easily met by choosing a finite volume. 
The domain of the mapping T is the Sobolev space of 
distributionally one-times spatially differentiable square- 
integrable potentials [10j[lTJ . This space is similar to the 
well-known domain of, e.g. the Laplace operator V 2 . As 
the domain is dense in the set the square-integrable func- 
tions we can uniquely extend our mapping via Cauchy- 
sequences to this set. In this way we can then also converge 
to Coulombic external potentials since they are square in- 
tegrable, similarly to the case of the Laplace operator. 
The space of the potentials [2] is isomorphic to the Ba- 
nach space of strongly measurable functions from [to , T] 
onto the Hilbcrt space of square-integrable functions in 
space with the norm ess sup to<t<T ||-y(i)|| 2 < oo, where 
ll u (£)l| 2 — fn^r v(rt) 2 and ess supp denotes the essential 
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suprcmum |10) . If the density n(rt) in ea.(|16[l is not diver- 
gent at any point then for all possible v under considera- 
tions n(rt)'S7v(rt) is square-integrable. If we then require 
the second time-derivative of the density d?n(rt) to be in 
the dual to the Sobolev space of potentials^, i.e. a require- 
ment which is slightly weaker than square-integrability, 
then from eq.© we see that the divergence of the local 
forces are also in this dual space. We thus restrict our- 
selves to finite densities having a second time-derivative 
in this dual space. Similar restrictions have to be posed in 
[5] as has been discussed in [2}UT] . With these restrictions 
the Sturm-Liouville problem (I16[) can be solved uniquely. 
The stronger restriction of square-integrable second time- 
derivatives of the densities and square-integrable internal 
force divergences is usually assumed to hold for all physical 
wave- functions [BJ. 

We start by deriving inequality eq. ([2"0"[) for the mapping 
V of Fig. ([1]). The main ingredient is the fundamental the- 
orem of calculus [15] . Suppressing the time-arguments we 
can write 

OH - OK] = J o dX^-[v + avi-v )]\ £=x . (22) 

for any operator expectation value 0([i>],i). In the case 
that we take = q(r) this equation yields 

q([vi],rt) - q([v ],rt) 

dt' f dr'x(rt,rY)( Vl (rY)- Vo (rY)) (23) 

to Jfl 

where we defined 

X (rt,r't') = -i f dX^ \[q Hx (rt),n Hx (r't')}\^ ) (24) 
Jo 

and Oh x is the operator O in the Heisenberg represen- 
tation with respect to Hamiltonian H\ with potential 
v\ = Vq + X(vi — vq). Equation (|24l) can be obtained di- 
rectly from the TDSE by evaluating the expectation value 
g([uf],rf) to first order in £ around A which amounts to 
linear response theory pQ. In the following we fix v\ and 
and derive an approximation for the induced difference 
in the local-force divergences. The linear response kernel 
X depends on this choice of potentials, i.e. \ = x[ v o>Vi\- 
Equation (123|) has the form 

f(rt) = f dt' f dr' X (rt,r't')g(r't') = (x<?)(rt). (25) 
Jto Jn 

One can then derive that 



||/(t)|| 2 <0(t) 2 / dt'\\g(t')\ 



(26) 



In this expression the function C(t) (the operator norm) 
is defined as 

g(*) 2 =sup l ! U9 . m L - (27) 



g ^o j; o dt>\\g(t>)\\ 



X A thorough discussion of these mathematical details can be 
found in [Tl] 



Note, that we have fixed x = xNi i] while ranging over 
different functions g. We assume C(t) to be finite, i.e. the 
operator norm of x to exist. Due to the dependence of x 
on the potentials chosen also C(t) depends on this choice. 
Using a Cauchy-Schwarz inequality the function C(t) can 
be shown to satisfy 



C(t) 2 < max / dvdr'x{vtyt' f. 

t'e[t ,t] J n 



(28) 



The integral on the right hand side of inequality (]26p can 
be manipulated as follows 



' II „/V\l|2 



dt' \\g(t 



it, 



^' e -«(*'-*o) e «(t'-*o)|| 9 ^/)||2 ( 2Q ) 
ft , , . e a(t-t ) 

< \\g\\l, / dt'e^ -*o) < M * * 



where a is an arbitrary positive number and where we 
defined 



\9\\l,t = ess sup (llff(f)lr' 

t'e[*o,t] v 



2 a (t'-t ) 



With this result inequality ([26]) can be written as 



II/IIL< 



C(tf 



2 

\a,f 



(30) 



(31) 



We will consider functions on an arbitrarily large but finite 
interval [t ,T] with T > t . If we define ||/|| Q = ||/|| q ,t 
and C = C(T) and apply this to eq. (|2"3"]) we obtain eq-ipD]). 
This concludes the first part of the proof. The discussion 
so far was related to the mapping V . Let us now discuss 
the second mapping V. 

We consider the Sturm-Liouville operator Q = —V - [nV]. 
By partial integration we find for two external potentials 
and v\ within the standard inner product 



(v \Qvx) - (Qvq\v\) = 



(32) 



dS ■ n(rt) (wi(rt)Vuo(rt) - w (ri)V«i(rt)) • 

on 

Hence, in order for Q to be self-adjoint the boundary term 
has to vanish. If n(rt) > e > in ft we know that we find 
solutions to the general inhomogeneous Sturm-Liouville 
problem [11] with the boundary conditions v — on dfl, 
making Q a self-adjoint operator. We then have an or- 
thonormal set of eigenfunctions [15 . If the density is zero 
at the boundary the same conclusion holds under the as- 
sumption of n^ 1 being square-integrable |14) . In general 
a boundedness condition at the edge is assumed to sin- 
gle out a unique set (a famous example is the Legendre 
equation which would correspond to an operator Q with 
density n(x) = (1 — x 2 ) on [—1,1] ). For an extensive 
discussion of these issues in the one-dimensional case see 
[16] , In general we have a set of orthonormal eigenfunc- 
tions {(/)i(rt)} with Q(pi — \i4>i and positive eigenvalues 
< A < Ai < since A; = (4>i\Q(f>i) = (V<fo|raV<fo) > 0. 
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The eigenfunction to the eigenvalue zero is cf>o = c(t). If 
we now subtract eq. (|16[) for V\ = ^[vq] from the one for 
v-2 = we obtain 

- V ■ [n(ri)V(«3(rt)-t>i(rt))] = g([«i],rt)-«([«o],rt). 

(33) 

We can then expand V2 — V\ and C, = q[v\] — q[i)o] in 
terms of the eigenfunctions of Q, i.e, (1)2 — vi)(vt) = 

X^o u *( i ) < M rt ) and C( r *) = X^o Ci (*)&(!"*). Since ( 
is a divergence J d 3 r((rt) = and therefore the constant 
function (j)Q does not contribute to the expansion. Like- 
wise the gauge fixing allows us to exclude d> from the 
expansion of vi — v\. Inserting the expansions into eq.(|33jl 
yields 

£•(*) 



Ui{t) 



Ai(t) 



(34) 



where Xi(t) > and with this we have 



Mt)-vi®\\ 



E 



1 LXJ 

v J i — 1 



A,-(t) 

q[vi](t)-q[v ](t)f. (35) 



If we now multiply eq. (|3"5)) with e~ alyt ~ to ^ and take a max- 
imum over the interval [t , T] we arrive at eq. (|2ip in which 
D = max te / {Ai(i)~ 2 }. This then together with eq. (|2"0|) 
establishes our main inequality (fT9)l . 

This equation can now be used to prove the uniqueness 
of a solution to eq. (|13p . Suppose we would have two fixed 
point solutions u and v, i.e. u = J-[u] and v = F[v\. 
Then by choosing <Ja — 2 CD for this pair of potentials 
in eq. (fT9l) we find 



\\v- 



= \\F\v\ - T\u\ 



< ^\\y- u \ 



(36) 



from which we conclude that \\v — u\\ a = and thus we 
have u = v. Hence, if a solution to eq.flT3]) exists then it is 
unique. This conclusion is equivalent to the Runge-Gross 
theorem 0]. It states that a density n([v],rt) cannot be 
produced by another potential u starting from the same 
initial state. This is now proven without assumptions on 
the Taylor-expandability in time of the potential. Suppose 
that the density of an interacting system is representable 
in a non-interacting system then this theorem guarantees 
that the effective potential producing the same density in 
this system is unique. This establishes the uniqueness of 
a Kohn-Sham scheme. 

Let us now address the existence of a solution to eq. ([TBf . 
This is a u-representability question for a given density. 
We see from ea.([57]l that the constant C = C(T) in eq. ([2"6|) 
is dependent on the response function \ an( i hence via 
cq. (|2"3"|) on potentials and V\, i.e. C — C[vq,v\\. We 
assume that a constant Csup = sup„ C[vq, F[vq\] exists 
when we range over all potentials vq. We now follow ex- 
actly the same reasoning as in the proof of the Banach 



fixed point theorem [TS]. Let — ^[vq] denote the k- 
fold application of the mapping J 7 on a given initial po- 
tential Vq and choose y/a > CsupD. Then cq. (fT9|) with 
a = CsupD/y/a implies \\v k+1 - v k \\ a < a k \\vi - v Q \\ a 
which means that the v k are a Cauchy series. Since the 
set of potentials is a Banach space with the norm (|30f 
[T01U3] and therefore complete this series converges to a 
unique v , i.e. v k — > v for k — > 00. According to our as- 
sumption the response function of eq. (|24|) exists and hence 
q is functionally differentiable and consequently continu- 
ous as a functional of v. Therefore lim/ £ _ ! . 00 q[v k ] — q[v] 
which means that v solves eq. (|13p and hence is a fixed 
point. This establishes the existence of a Kohn-Sham sys- 
tem corresponding to the density n in ea. (|13[) provided 
there is a supremum sup„ C[vq, ^[vq]] when we range over 
potentials vq in a non-interacting system. 

We have generalized the uniqueness and existence the- 
orems of TDDFT to cover the case of external potentials 
that are not Taylor-expandable in time using a new formu- 
lation of the problem as a fixed point problem in terms of a 
suitably chosen norm. We showed that a potential gener- 
ating a given density is unique provided that the response 
function of eq. (|24[) has a finite operator norm and the 
second time-derivative of the density to obeys a slightly 
weaker condition than square-integrability. We further 
showed the existence of an external potential producing 
a given density provided there exists an upper bound for 
the operator norm of all response functions between suc- 
cessive potentials of the iterative procedure. The main 
idea has several applications in other research fields. For 
instance, we will extend the ideas presented here to other 
density functional theories, e.g., current-density functional 
theory [17] and also study the convergence properties of 
optimized effective potential approaches [T8H20] . As the 
presented iteration scheme poses a way to deduce an ex- 
ternal potential generating a given density, this could also 
prove useful for quantum optimal control theory [21] . 
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